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Abstract 

The quark-gluon mixed condensate g^cr^F^q) is calculated in 
the Gaussian approximation of the Field Correlator Method. In the 
large N c limit and for zero mass quarks one obtains a simple result, 
m o = 9 ^ a ^^ l ' q ^ = where a is the string tension. For a standard 
value a = 0.18 GeV 2 one obtains uiq = 1 GeV 2 in good agreement 
with the QCD sum rules estimate m 2 , = (0.8 ± 0.2) GeV 2 and the 
latest lattice result tuq = 1 GeV 2 . 



1 Introduction 

The mixed quark-gluon condensate (QGC) is an important characteristics of 
the nonperturbative QCD vacuum, which together with the quark condensate 
(qq) signals the chiral symmetry breaking. Moreover, the QGC measures the 
average interaction of the quark color-magnetic moment with the vacuum 
fields, which is an important ingredient of the quark dynamics in the vacuum 
(e.g. it is this term which gives attraction of in quark zero modes). 

In the QCD sum rules the QGC plays an important role PQ and the 
phenomenological analysis suggests the value of m 2 , in the range = (0.8 ± 
0.2) GeV 2 pQ, see j2j for a review. One should stress at this point that 
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for a nonzero quark mass m the (diverging) perturbative part should be 
subtracted. 

As will be seen below the resulting nonperturbative dependence of mg 
on m is very weak in agreement with lattice data. Lattice studies of QGC 
[H] , [I] , have not yet converged to a definite prediction. A problem there 
is the extrapolation to zero quark mass and the quenched approximation. 
In ref. [1] the simulations are done in the quenched approximation, the 
condensate is measured by use of staggered quarks, and the result for Wq is 
definitely larger than the sum-rules value. Ref. [5] uses an optimized version 
jH] of domain wall fermions, which are better in principle for the chiral limit, 
again in the quenched approximation. Their result is m\ = lGeV 2 , which 
agrees with QCD sum rules. It is therefore worthwhile to calculate QGC by 
a different nonperturbative method. 

In the framework of the Field Correlator Method (FCM) [7] the color- 
magnetic quark gluon interaction term ga^F^ enters essentially in the Fock- 
Feynman-Schwinger Representation (FFSR) of the quark propagator in the 
vacuum background field In particular the quadratic average of this 
term defines the hyperfine qq interaction where the nonperturbative part is 
proportional to the field correlator {F lu/ (x)F pa -(0)) measured on the lattice 0. 
Even more important the term ga^F^ is in the contribution to the bound 
quark self-energy [121, where it is of paramagnetic character, i.e. negative 
and strongly decreases the masses of hadrons, putting them in accordance 
with experimental data ^T]. Explicit correction to the bound quark mass 
squared is (TU] 

Am = rj (1) 

7T 

where r\ = rj(mT g ) is a calculable function of the quark current mass m, 
renormalized at the scale of 1 GeV. The function 77 is given in ^U] and in 
the Appendix below and for zero quark mass is normalized to one: f](0) = 1. 
We calculate in the next section the QGC, or rather the parameter m\ in the 
same way, as it was done in jTU] for Am 2 , with the result 

9 A . 9 16(7 . . 

mo = ~4Am 2 = 77. (2) 

For a = 0.18 GeV 2 one obtains ttiq = 1 GeV 2 which is in agreement with the 
lattice data 0, and with the QCD sum rules estimate quoted above. 
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2 Calculation of ml 

We proceed in the Euclidean space-time and write 

(qg^F^q) g ,A = tr(ga IML ,F fML ,(x)S g (x, x)) A = tr(S q (x, x)ga fMU F fMU (x)) (3) 

where S q (x,y) is the Euclidean quark propagator, for which one can write 
using the FFSR 

S q (x, y) = (m + D)~l = (m - D) x (m 2 - D 2 )~l = 

poo rs 

(m-D) ds(Dz x>y e- K $ z (x,y)P F ex V X{z{r))dr. (4) 
Jo Jo 

In (JIJ) the following notations are used: K = m 2 s + | J S i^dr, = d p — 

igA^, (Dz) Xj y is the path-integral measure for paths starting at y and end- 

ing at the point x, {Dz) x , y = lim^ U^ =1 (^$^) $fce*<£^x M-M-iM 

while is the phase factor (parallel transporter) along the path z^{r) 

®z(x,y) = -P/iexpzg Jy A^dz^, with Pa,Pf - the ordering operators of the 

matrices and A(z), where X(z) is defined to be 1 

K z ) = 9<Ti*>Ff*>{z)\ °V = jt(7m7* - 7*7/*) ■ ( 5 ) 

For what follows it will be advantageous to take in © A(z(r)) = g^a^F^z^r)) 
since the functional derivative t-t-t at r — ► or r — ► s inside the FFSR 
(J3J) brings down additional factor A(y) or A(x). When one has y = x, as 
in 0, then both contributions add, which formally is obtained by putting 
g(0) = g(s). In this way one can rewrite (J3J) as follows 

(q(x)X(x)q(x)) = tr(X(x)S q (x, x)) = 2trj-—r(S q (x, x)). (6) 

As the next step one can write the average (S q (x,x)) in the form of cluster 
expansion jjj 

(S g (x,x)) = ((m-ip) J dse~ K (Dz) xx exp{~ J dv Xp J dv au (gF Xp gF au )} 

°__ (7) 

1 the definition of in (5) (as well as in |7],|H]) differs from the standard definition 
in QCD sum rules, where enters | instead of \ in (5). Therefore one obtains additional 
factor 2 in the definition of m\ in (17). 
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where only the contribution of the lowest cumulant (FF) is retained in ac- 
cordance with estimates [T2] . and the nonabelian Stokes theorem is used to 
express A p through F pu , with the notation 

dv Xp = ds X p - icrxpdr, gF Xp dv Xp = 

= gF X p(u)ds X p(u) - ig(T)a X pF X p(z(r)) (8) 

and ds Xp is the element of the area of the surface enclosed by the contour 
z m (t), z M (0) = Zn(s) = Xp. Performing differentiation in © one gets 

(qXq) = 2g 2 <jp U cr X p ds(Dz) xx e~ K (m - ip) dT(F Xp (u(T))F p „(x))x 

x exp{-y J dv Xp J dv (TV {F Xp F (TP )}. (9) 
Using the identities |8| 

roc rs roc roc 

(Dz) xx = (Dz) xu d%(Dz) ux , ds drf(s,r)= ds dr/(s + r,r) 

Jo Jo Jo Jo 

(10) 

where f(s, r) is an arbitrary function, one has 

(qXq) = 2a pu a Xp J (G(x, u)S q (u, x)) D x 2) p41l/ {u - x)d A (u - x). (11) 
Here we have defined as in [7j 



D \p,pu( z ) = {b\pbpv-b\vbpp)D{z)+-(d x Zp5p U +dpZ u 5 X p-d x z u 5pp-d p Zp5 Xu )Di(z) 

(12) 

and 

G(x,u) = J dre~ K (Dz) xu exp{-^ J dv Xp J dv au (gF Xp gF au )}. (13) 

Note that Gq(x, u) and S q (u,x) share common factors depending on a piece 
of common z between u p and x p and in general do not factorize. 

At this point we shall use the properties of the correlators D(z), Di(z) 
found on the lattice 0, in the quenched case one has 

D{z) = 3£>i(z) = D(0) exp(-\z\S), 5 = 1/T g « 1 GeV. (14) 
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Analytic calculations based on the gluelump spectrum [12] suggest even 
larger value, 5 ~ 1.4-=-1.5 GeV. The string tension a can be expressed through 
D(z) (the correction due to higher correlators is limited by the Casimir scal- 
ing arguments to a few percent [T2]) 



a 



\ J D(z)d 2 z. (15) 



Since the distance \u — x\ is of the order of T g , we can now use the argument 
of the small T g limit (large 5) for the constant a to factorize the product 
(G(x,u)S q (u,x)) as follows 

\im^(G(x, u)S g (u, x)) = Go(x — u)(S q (x, x)). (16) 

This approximation is equivalent to the expansion in the parameter £ = 
crT 2 <C 1. As the result one obtains the following representation for the ratio 



ml = 2 



(qXq) 



^ v a Xp J G (z)D^(z)d 4 z (17) 



° " " (qq) 

Go(z) is easily calculated using (JT3|) to be the free propagator of the scalar 
quark with mass m, 

Go(z) = n*Hf!) (18) 



where K\ is the McDonald function, and m is the current (pole) quark mass 
normalized at 1 GeV. 

Taking into account that 2 



one obtains for wiq 



poo 

m 2 Q = 12m / z 2 dzK 1 (mz)(D(z) + DAz)) (20) 
Jo 

or, with the help of (JHj), 

/•OO 16<T 

16m / z 2 dzK 1 (mz)D(z) = <f(m/S) (21) 

JO 71 



ml ~ 



2 note the misprint in Eq.(15) of [9], where coefficients of D,Di differ from those in 
(|19J) . Nevertheless the final result in Eq.(29) of ^01 is the same as in our Eq. Q due 
to the relation D\ m -|D [Hj valid for the quenched case, considered here, whereas in the 
unquenched case one obtains instead of (1): Am^(m — * 0) = —3 J °° zdz(D + Di) = — ^er. 
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where we have defined 



/ 

Jo 



oo 



(p(m/5) = m5 2 



z 2 dzKi{mz) exp(— 5z); <£>(0) = 1. 



(22) 



It is easy to see with the help of (j!5|) that in the limit of small quark 
mass, m — > 0, one obtains for a = 0.18 GeV 2 (in the quenched case) 



It is appropriate at this point to discuss the accuracy of our result (|23(l. 
The main uncertainty appears in expressions ()14j). (jTSjl and (fTB|) and we 
consider the accuracy of the corresponding approximations point by point. 

The lattice calculations [9] of D(z) and D\(z) define the amplitudes A, Ai 
and slopes 5,5\, the first ones are reabsorbed in the value of a, while the 
latter are equal with accuracy of few (1-2) percent to the value given in 
()14)l . The approximation of ()15jl reduces to the neglect of higher correlators, 
contributing to the observed string tension a. This accuracy was tested in [12] 
using the Casimir scaling and is also of the order of few percent. The largest 
possible error may come from the replacement (jTEJ), where one can use the 
fact that the integral over d 4 (u— x) in (fTTj) is taken with the weight D&'iu—x). 
The latter is exponentially decreasing at the distance 1/5, while the range 
of G(x,u) is defined by the confining exponent in (JT3J), which produces the 
effective quark mass, computed through a and equal to 0.35 GeV for the 
lowest state (see [11] for references and explicit calculations). Introducing 
this mass instead of m in (JTBJ), (|2Uj) . (|2*T|) one obtains <p ~ 0.75-=-0.8, and using 
(|2*T|) one comes to the conclusion that mf) is in the range 0.7 GeV 2 < m 2 , < 1 
GeV 2 . This range lies very close to the limits predicted in the QCD sum 
rules. 

The explicit analytic form of <p(x) was obtained in jTU] and is given here 
in Appendix. For 5 = 1 GeV, and m = 0.175 GeV, 1.7 GeV and 5 GeV one 
obtains respectively tp = 0.88, 0.234 and 0.052. 

The resulting value of m 2 , (J23)) is in agreement with the QCD sum rule 
estimates j2], and with the lattice evaluation of tUq, namely tUq ~ 1 GeV 2 in 
0. One should note, that there is a large perturbative contribution to m 2 , 
for nonzero quark mass m proportional to mAfjy ~ m/a 2 , which should be 
subtracted to get agreement with purely nonperturbative result (123(1 . 

On the other hand the purely nonperturbative behaviour of TTZq ctS 8l func- 
tion of the quark mass m, or rather the ratio t = m/5 is given in appendix, 



mg(m — > 0) 



—o = 0.92 GeV 2 . 



(23) 



7T 
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Eq. (A.6), 

mlit) = — (1 + t 2 (4 - 3 In \) + 0(t A )). (24) 

7T t 

The values obtained from (J24|) agree well with the lattice measured 

values in [3] for ma > 0. Indeed for three values of ma, ma = 0.05; 0, 1 and 
0.15 one obtains from (J2U) taking a = 0.18 GeV 2 , and a" 1 = 1, 979 GeV 0, 
wiq = 0.434; 0.393 and 0.342 GeV 2 respectively. This should be compared 
with the values m§(ma) measured in |3J and equal to 0.371; 0.311 and 0.290 
GeV 2 . At the same time the limiting extrapolated value m\{ma = 0) ~ 1 
GeV 2 obtained in f5], agrees with the theoretical one, given by Eq. (J2U), 
m\(ma = 0, theory) = 1 GeV 2 . One should have in mind, that chiral quark 
mass corrections present in both the quark condensate and the QGC are 
cancelled in the ratio m\ to the leading order in cT 2 , so the remnant ma 
dependence in m\ comes from quadratic terms in (J24|) and linear perturbative 
terms mentioned above. 

Recently a study of thermal dependence of m\{T) has been reported in 
|14j . where m 2 , was found almost independent of T up to T = T c . This is in 
general agreement with our expression (23), since a is roughly constant in 
that region, but more detailed check of behaviour near T c is desirable. 

Summarizing, we have obtained a simple nonperturbative estimate for 
the ratio of condensates, which is in a reasonable agreement with the QCD 
sum rule results, and lattice results in [3] for nonzero ma and zero ma limit. 
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collaboration at the first stage of this work. 
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Appendix 



The function <p(t),t = m/5, defined in Eq.(22) can be written as (note 
the difference in definition here and in [TU] ) 

poo 

<p(t) = t / z 2 dzK 1 (tz)e~ z (A.l) 
Jo 

where K\ is the McDonald function, K\[x){x — ► 0) ~ -, so that for t — 
one obtains 

<p(0) = 1. (A.2) 
For t > the integration in (jA.l|) yields two different forms; e.g. for t < 1 , 

3t 2 , 1 + y/T^¥ l + 2t 2 . , 

= - (TT^ ln i + (73^ 

while for t > 1 one has instead, 

= - (^2^1)5/2 arctan(v / ^T) + (A.4) 
For large t one has the following limiting behaviour, 



For small t one obtains expanding the r.h.s. of (|A.3|) 

pit) = 1 + t\A - 3 In - t ) + t\ 7 - - y In |) + 0(t 6 ). (A.6) 
Some numerical values are useful in applications. 

</?(0.175) £ 0.88, p(1.7) = 0.234, y>(5) S 0.052 
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